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SUMMARY 


Grumman, under contract NASl-14643, developed the computer code CAVES 
("A General Transient Heat Transfer Computer Code Utilizing Eigenvectors and 
Eigenvalues") that performed a specialized transient thermal analysis. Under the 
present funding a new extension of CAVES was developed which automates the creation 
of the input. 

The new software- 

• Utilizes Tektronix 4014 Graphic Scopes to display models 

• Utilizes Tektronix Tablet Digitizers to generate model geometry 

• Graphically displays the model geometry 

• Creates a finite difference digital dataset with all conduction links 
computed ready for input to CAVES. 
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Section 1 


INTRODUCTION 

The Network Generator software created for Langley Research Center automates 
the creation of finite difference thermal models for input to the CAVES code provided 
Langley under contract NASl-14643. The software makes use of the interactive 
capabilities of the INTERCOM CDC timesharing system on the Langley Research 
Center CDC 6600 computer. It utilizes the Tektronix 4014 Graphic Scope and digitizing 
interface and makes use of the appropriate Tektronix software resident in the Plot-10 
package. The Network Generator utilizes the graphics hardware to expedite the 
transformation of geometric data to digital data and provides visual confirmation of 
geometric and nodal networks. 

The Network Generator is broken into two distinct programs : geometric data 
mput and model generation. The geometry data mput software is labeled NNDIG and 
makes use of Tektronix digitizmg hardware to create a data set containing the pertment 
geometric mformation needed. This option provides the user with the capability to 
transfer geometric data from a scaled diagram displaying the model to be analyzed to 
the computer as a digital data set. This software does not have to be run if the 
user wants to generate the resultant geometry tile by hand (if digitizmg hardware were 
not available) . 

In either case, the geometry is fed to the model generation program (NNGEN). 

It IS here that the actual fmite difference model is created. NNGEN uses the geometry 
supplied to it with the specific properties of the materials to be analyzed and creates 
the digital data set that represents the finite difference model. Software establishes the 
numbering sequence, determmes all geometric voids withm the model, and bookkeeps 
all couplmgs. The software can generate two- or three-dimensional models and 
provides the user with the ability to have mternal or boundary (surface) nodes. 
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Mr. James L. Hunt of the High Speed Aerodynamics Division, Langley Research 
Center, Virginia, served as the NASA technical monitor for the program. 

At Grumman, the contract was admmistered by the Advanced Development 
office under Mr. Fred Berger, Manager of Advanced Development System Engmeermg, 
Dr. Keimeth A. Rathjen was Study Manager and Mr. Joseph V. Palmier! was 
Engineering Specialist, Software Development. 



Section 2 


NODAL NETWORK GENERATOR DESCRIPTION 

The Nodal Network Generator package was created to automate the procedure of 
creatmg a digital math nodal to be used as input to the CAVES code (Reference 1). 

The package consists of two distmct programs; a Nodal Network Digitizmg (NNDIG) 
program and a Nodal Network Generator (NNGEN) program. The package makes use 
of the Tektronix digitizing tablet and graphic scope (Model 4014) to transfer geometric 
nodal information directly to the computer and plot the resulting nodal network. 

The program enables a user to create a math model by merely digitizing the 
vertices of homogeneous sectors and digitizing a grid increment that defines the nodal 
surface. The program will automatically compute the geometric centroid of the nodal 
areas and place the nodes there , then automatically compute all the nodes each node 
sees , the conductance between nodes , the capacitance of the node , and the radiation 
area of the node surfaces in the x, y, and z direction. The program will determine 
whether the body has voids within it and will establish boundary nodes to identify 
this condition. 

The program has two nodal options that can be employed; nodes can be placed 
at the centroids for all nodes or placed at the surface for boundary nodes with the 
conductances adjusted for the distance from the surface to the centroid. The program 
also provides a three-dimensional feature which allows the user to create a three- 
dimensional regular model. The sections to follow will explain all of these features 
in detail. 

Once the software has been accessed and executed, the Nodal Network Generator 
will create an mput file capable of be mg run with CAVES. The only change to the data 
that may have to be made would be to mclude the boundary conditions using the areas 
created by the program, mclude interface conductance effects if desired, and in- 
corporate any radiation couplmgs. 
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Section 3 


PROBLEM SETUP & PROGRAM LIMITATIONS 

The first step in the development of a thermal math model is the creation of 
the nodal model on paper. Here the user can superimpose a nodal grid on the 
physical boundaries of the object to be analyzed. The Nodal Network Generator needs 
this first step as the engineer did in establishing a digital math model. The graphic 
representation of the nodal network must meet certain criterion before iL can be used in 
conjunction with NNDIG. Consider Figure 1, a structure made up of a number of 
different materials. Each block shown (I - VI) represents a sector. This will be 
explamed later. The geometry should first be laid out m a convenient orientation 
relative to a righthand Cartesian coordinate system. Because the Nodal Network 
Generator only analyzes regular three-dimensional shapes, we need only be concerned 
with the top view of the structure . 

The external boundaries of the object should then be drawn to scale (Figure 2). 
Once the perimeter has been established, the body itself should be broken down into 
sectors. Sectors are defined as homogeneous divisions of the body being analyzed and 
having isotropic properties withm the sector bounds. The sector can have no more 
than 15 vertices and all angles must be convex, i.e. , no internal angle can exceed 
180°. The body itself can have no more than 20 sectors. The geometry displayed m 
Figure 2 was broken up mto six sectors. All sectors have uniform thickness normal 
to the plane of the paper. Sectors I and II are made up of the same material (steel) and 
have identical properties as are sectors III and IV (aluminum) . Sectors V (steel) and 
VI (alummum) have unique properties of their own. Notice that although Sectors I and 
II or III and IV have identical properties and would have met the homogeneous and iso- 
tropic rule as one sector, they had to be broken up into two independent sectors 
because an internal angle exceeded 180°. 

With the sectors and sector vertices identified, it is now necessary to establish 
a network grid. The Nodal Network Generator is designed to set up rectangular and/or 
triangular nodes. The node boundaries are established by the grid size and its mter- 
action with the physical boundaries of the model. The network grid is a system of lines 
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FIG 1. STRUCTURE TO BE ANALYZED 

parallel to the x and y axes that delineates the node boundaries of the body. The grid 
lines must be established running parallel to the x and y axis and will ultimately define 
the actual node size. 

To establish a grid, the user must assure that: 

• Each vertex will be intercepted by a grid line in both the x and y direction 

• Each sector surface that is intercepted by a grid Ime must also be 
intercepted at that point by the orthogonal grid line. (This ensures 
that sloped boundaries of a sector are divided into triangular nodes 
and not trapezoidal nodes) 
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I I 

(0,0) — (30,0) 

2165-002B 

FIG. 2 MODEL LAYOUT 


• When employing the boundary option (Section 5) , each sector must be 
divided at least once by a grid Ime in the primary plane of the body 
(x-y plane in Figure 2). This allows the program to '^decide*^ where 
the boundary is located and, thus, where to place boundary nodes. It is 
mconsistent with the boundary option to allow a node to "see’’ boundary 
condition on opposite sides of a node. (This does not apply to the 
direction normal to the primary plane of the body because a 
one-layer model can be used with either option). 
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Figure 3 displays the step-by-step development of a nodal network. Note the 
incremental steps defined by the grid line locations in the x and y direction. All 
grid Imes shown that have an asterisk indicate a grid line position that must be m 
the model. All non-asterisk lines are optional and at discretion of the user. Because 
it was stated previously that each vertex of a sector must have a grid Ime pass through 
it in both the x and y direction, when establishing the grid the user should first 
project Imes parallel to both axes emanating from each vertex as is demonstrated 
in Figure 3 by the arrows shown. Grid lines XI, X3 , X4, X5, X6, X9, Y4, Y6 , 

Y7, Y8, and Y9 were created this way. Note the Y7 projection crosses a skewed 
line of sectors V and VI, and, according to the rules previously stated, an orthogonal 
grid Ime must emanate from the intersection. This therefore established grid Ime 
X8. Lmes Y2 , Y3, and Y5 were put in to better defme the analysis by the user. 

Note that the user included X2 to satisfy criterion from boundary option. As shown, 
depending on the shape of the model analyzed, the number of nodes created may be 
large just from required grid lines. 

Once it has been established, the grid work (Figure 4) will represent the nodal 
boundaries of the model. Note by usmg this method all nodes are m line with one 
another, which is a necessary condition for the Nodel Network Generator. 

Now that the model is established on paper, the user can access NNDIG and 
NNGEN to create an mput data set to CAVE3 that contams all the network couplmgs 
and values of capacitance and conductance . 

The Nodal Network software will provide the user with the option of placing 
nodes at the geometric centroid of the node or at a predescribed surface location. The 
program also provides the user with the ability to create a three-dimensional model, 
in which case the network developed will be mirrored down to a select number of layers 
(up to 10 layers). The nodes will be mcremented and three-dimensional couplings 
will be established. All of these options will be described, but lead to another 
limitation. The program can currently create a model containing up to 500 nodes. 

Note the limitation on the amount of nodes created, maximum number of sectors, 
and maximum amount of vertices per sector are totally arbitrary and were limited 
only to reduce the amount of computer core needed to run the program. Any one or 
all of these limitations can be mcreased by simply re-dimensionmg the program. 
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Section 4 


INPUT DATA - MANUAL OR DIGITIZED (NNDIG) SETUP 

The input needed to create the digital thermal model is basically the descrip- 
tion of the physical geometry, i.e. the coordinates (x and y) of all sector vertices, 
and the x and y position of the grid lines and the material properties, conductivity, 
the specific heat, the density, and thickness of the material. The units inputted 
should be consistent, i.e. , if the scaled drawing is measured in feet, make sure the 
thickness is also m feet and the conductivity contams the units of feet, etc. The 
program contams a feature to automatically create this data set us mg the Tektronix 
graphic scope and digitizmg board. Although it is easy to use, this feature is not 
necessary. The user can manually create the data set needed as input. 

Figure 4 displays a nodal breakdown of the network to be analyzed. Figure 5 
displays the actual format of the data needed. Figure 6 represents a data set 
created manually. Note the program requires the number of sectors created within 
the model, m this case six (Figure 6). It also requires the number of layers and 
overall thickness of the model, which we input as 3 and 15. , respectively, i.e. , 
we chose to obtain three nodal cuts through the model and the overall depth of 
the model is 15 units. The next card is the respective thiclmess of each of the cuts, 

I.e. , 3. , 5. , and 7 units. If the cuts have identical thicknesses, the user could 
have inputted a zero (0.) in which case the program would automatically compute the 
thickness by dividing the overall thickness by the number of layers. As shown, three 
unique values of thickness were supplied. 

The program then needs the number of vertices tor the tirst sector (5) , the 
density (489), the specific heat (.1), and the conductivity (27). The next card should 
have the x and y coordinates of the first vertex, the next card the coordinates ot the 
second vertex, etc, until the user has mput the coordmates of each vertex on 
individual cards. It should be noted that the coordmates of the vertices must be entered 
m a contmuous clockwise pattern. Once the vertices are entered for sector I, then 
the property card is repeated for sector II. The user must supply the number ot 
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INPUT FORMAT 

16A4 

215, F104 
10F7 2 

12, 3F104 
2F123 

(IVERT+ 1) sector card -XVERT (I, IVERT), YVERT 

The next group of card defines the coordinates of the grid lines 
X(1),0 

X(K),0 
0 ,Y(1) 


1st Card — TITLE {64 characters) 

2nd Card - ISECT, NLAY, THICK1 

3rd Card- THICK (1), THICK (2), THICK (10) 

The next group of cards is required for each sector (I = 1 to ISECT) 

1st sector card - IVERT (I), DENS (I), CP(t), AK{I) 

The next group of sector cards define the vertices (J) for the sector (I) 
2nd sector card - XVERT (I, J), YVERT (I, J) 


0 ,Y(L) 


ISECT 

NLAY 

THICK1 

THICK (K) 

IVERT (I) 

DENS(I) 

CPU) 

AK (I) 

XVERT (I, J) 
YVERT (I, J) 
X (K) 

Y (K) 

2165 006B 


DEFINITION OF VARIABLES 

number of sectors 

number of layers 

overall thickness of model 

thickness of each layer where K = (1 NLAY) 

number of vertices in sector I 

density of sector 1 

specific heat of sector I 

conductivity of sector I 

X vertex number J of sector I 

Y vertex number J of sector I 

X coordinate of grid line where K = 1 to last X grid line 

Y coordinate of grid line where K = 1 to last Y grid line 


FIG 5. INPUT DATA FORMAT & VARIABLE DESCRIPTIONS 
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vertices contained in sector II, its density, specific heat, and conductivity. This 
card is followed by the coordinates of the vertices of sector III, and so on until you 
have inputted the material card and vertices for each sector of the model. Note when 
supplying the vertices , the user should repeat the first vertex of each sector as the 
last to close the body for plotting. 

The input needed now is the x location of the first grid line followed by a zero. 
This is continued until all x grid locations are defined. The same procedure is used 
for the y grid location, except that 0. is input first followed by the y coordinate of the 
grid location. The data file (Figure 6) is now ready for execution with NNGEN. 

The above data file could have been automatically created in a fraction of the 
time using the digitizing option of the program. Figure 4 will again be used but this 
time in conjunction with the digitizer. An axis had to be drawn on the figure with ref- 
erence point clearly denoted,^ as shown. Notice the grid tick points are marked on the 
respective axis. In this case we have nine grid marks in each direction. Note the 
number of grid points in the x direction does not have to equal the number of grid points 
in the y direction. The user should then affix the drawmg to the digitizing table and 
access the program by typing at the Tektronix terminal: NNDIG (Figure 7). 


nndig 

PLEASE ENTER A TITLE FOR THE DATA FILE (60 CHAR - MAX) 
3“d model for report 

ENTER THREEC3) REF. PTS. (XRF(I),YRF(I ),XRF(I ) ) 

AS DECIMAL NUMBERS SEP. BY COMMA 
POINT 1>0.,E6.,0., 

POINT a>0.,0.,0., 

POINT 3>30.,0.,0., 

DIGITIZE THESE THREE (3) REF. PTS. 


DIGITIZE TUO(E) PTS. (LOUER-LEFT & UPPER RIGHT) 
EACH HALF INCH AUAY FROM PICTURE FRAME 

2165-046B 


FIG. 7 SAMPLE TERMINAL SESSION NNDIG 
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The program then asks the user a number of prompts, the first being, "Enter 
a title for use as the first card in the data file," The program then calls for the 
reference points. As shown, the three points marked on Figure 4 were input at the 
terminal and then were digitized in the order they were input. The user merely placed 
the digitizing cursor to each point and recorded the data. He then placed the cursor 
on the remaining points and did the same. Once the three points were recorded, he 
pressed the carriage return to transmit them to the computer. This establishes the 
reference of the drawing, and any point recorded within these reference points will 
have a value based on these points. The user is then requested to digitize (same 
procedure as above) a pomt to the lower lefthand comer of the drawing outside the 
reference frame and the upper righthand corner. This is used to scale the drawing to 
the screen size and sample marks are denoted in Figure 4. Once accomplished the 
screen will be erased. 

The program will then ask the user how many sectors are to be digitized, how 
many layers the model will have , and what is the overall thickness (depth in the plane 
of the paper) of the model. As shown in Figure 8 the model has six sectors and will 
be a three-dimensional model consisting of three distinct cuts with an overall thick- 
ness of 15 units. The program then asks the respective thickness of each layer. As 
shown, the layers will have unit thicknesses of 3. , 5. , and 7. , respectively. If they 
all had the same thickness or there were only one layer, a carriage return would have 
been sufficient and the program would automatically establish the correct thicknesses. 
This data is now used by the program to create the second and third cards of the data 
file. 


Diror TIC NUMER OF NonoosNOis senKt TO m oioiriso 

AND lUaCR OF lAVEIB MO IWGRnLL THICKNESS 
>€.3.15.. 

iENTCR THICKNESS OF EACH LATER (liFT.O) 

C.R. ALOC INDICATES UHlOAn THICKNESS 
>3.,S..7.. 

2165-008B 


FIG. 8 INPUT OF OVERALL PARAMETERS 


17 




The program then queries the user (Figure 9) for the density, specific heat, 
and conductivity of the sector to be digitized. Once recorded, the user can begin 
digitizmg the vertices of the sector. The vertices must be recorded in a clockwise 
manner starting at any vertex. The user should (for plotting purposes) return to the 
vertex he started with to close the geometry, i.e. , the first vertex will also be the last 
vertex recorded. Once it is digitized, the geometry should be transmitted (carriage 
return) . If there are more than nine vertices , you should transmit (carriage return) 
at or before nine , then proceed to record the remainder , never exceeding groups of 
nine before transmitting. Once accomplished, the user should press a backspace once, 
the carriage return to indicate he is finished with this sector. This procedure is 
repeated until all sectors are digitized (Figures 9 through 14) . 


kwrcf OCNBITV. CP. CNO eOCUeTIVITY OF KCTOI 1 
loMITtS HC CBITtCCS OP SECTOR I 



FIG. 9 DIGITIZE SECTOR I 
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The program will then ask the user to digitize the lines defining the subdivision 
along the x axis (Figure 15) that have been laid out on Figure 4. To do so as before, 
each point is recorded usmg the digitizing cursor and transmitted with the carriage 
return. Groups of up to nine recorded points may be stored before transmission is 
required (carriage return). Once complete, the user should press a backspace and 
carriage return to indicate this portion is finished. Figure 16 displays the same 
procedure for the y increments. Once completed (backspace then carriage return 
pressed) , the screen will erase and a plot will be drawn of the sectors digitized 
(Figure 17). When the user is finished looking at the plot, a carriage return will 
erase the screen, and a data set will be established as shown m Figure 18. 
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ENIE^ OCNSm. CP. MO eONDUCTIOITY OT KCtOR 4 
>i79.,.2a.tat.. 


oraiTis TIC uorriccs of sector 4 
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FIG. 12 DIGITIZE SECTOR IV 


ENfl^ OOOITT. CR. MO CQMUCTIMITY OF XCIOI S 

>4n.,.t,27.. 

0I6IT1Z TIC OERTICEB OF SECTOR S 

i 

t 


2165-0138 

FIG. 13 DIGITIZE SECTOR V 
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EMI^ OEMSITY. MO CCMOUCTIUITY OF SBTDll I 
>41).,. 1.3?., 


OI«ITI2E TK lOnCES OF 9ECTDR • 
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F G. 14 DIGIT ZE SECTOR VI 



FIG. 15 DIGITIZE X GRID 







>MM 0tCtTI2E H« NQML tNOODITB » >Y- 
sfMinNO mm niN -r- to mu incubino bounmib 
incss OMKiMc fcniM wem remt to mam data : 

> 
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FIG. 18 SAMPLE INPUT TO NNGEN (CREATED BY NNDIG) 
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Section 5 


MODEL GENERATOR (NNGEN) 

Once it has been created, the geometry data file must be fed to the model 
generation portion of the program where the user has two options he might employ. 
The user can create a nodal network with the node centers located at the geometric 
centroid or the user can opt to have all boundary nodes moved to the surface in a 
predetermined fashion that will be described later. Each option will culminate in the 
creation of a data set capable of runnii^ with the CAVE 3 code. 

Using the geometry data file just created (Section 4) , the user can create a 
thermal model by accessing the network generation program (NNGEN). As demon- 
strated in Figure 19 , the program was accessed and it queried for an input of 
tolerance , a mechanism provided for in the software to allow for inaccuracies when 
digitizing. 


nngen 

INPUT TOLERANCE E1E.5 

♦ 1 

ENTER NODE N0., COND N0., RAD N0«, YOU UISH TO START UITH. IF 1 CR 

DO YOU UISH TO SEE A PLOT OF THE INTERNAL NODE CONFIGURATION? 

yes 

2165-018B 

FIG. 19 SAMPLE TERMINAL SESSION NNGEN 

The routme establishes its bookkeepmg scheme by looking for node centers that 
are m Ime with one another, i.e. , common x or y coordinates. When digitizing, it is 
extremely difficult (and most unlikely) for the user to repeat a point exactly. An 
example is shown m Figure 18. The data set represented here is one that was pro- 
duced by the digitizmg portion of the program. Note any set of vertices of any sector. 
Most sectors will show repeat coordinates, i.e. , a boundary that is parallel to either 




the X or y axis would create like x or y coordinates in two vertices. Note the first 
sector shows like x coordinates of 3.967 and 3,954; close but not exact. This error 
(0,013) was caused by the inability of the user to repeat exactly the x coordinate using 
the digitizing cursor. Although it is small, the error would have misaligned node 
centers causing the bookkeeping scheme to fail. 

To overcome this , the program requires an input of tolerance which is used 
as a ± buffer on node center to determine allowable node alignment. The tolerance 
has been made an input because it depends upon the scale of the geometry and node 
size. In any case, it should always be less than half the distance of the smallest x 
or y nodal increment. Because the tolerance is only used to bookkeep the orientation 
of the nodes , it does not alter the actual value of the node centers and therefore does 
not vary the dimensional parameters. 

A value of 0.1 was inputted for tolerance (any value up to approximately 1.0 
could have been used in this case) . The program then questions the user as to a 
starting node number, conductance number, and radiation number. This option was 
provided to allow the user to build a model in steps. In the terminal session provided, 
the user entered a carriage return indicating he wanted all numbering to start at 1. 

The routine then establishes the nodal arrangement of the geometry . Starting from the 
lower lefthand comer of the drawing, the routine establishes consecutive numbering 
of nodes and identifies any voids in the geometry, as demonstrated in Figure 20. The 
rectangular node centers are placed at the geometric centroid of the nodal area. In 
the case of triangular nodes , the node center is temporarily placed at the midpoint 
of the hypotenuse. This is only done for bookkeeping purposes for node alignment. 

At this point the software has established the bookkeeping of the model, i.e. , it knows 
which nodes see which nodes, which nodes see void (boundary) for all three dimensions. 

The program then asks the user whether he wants to see a plot of the internal 
node configuration, i.e. , with all node centers located at the geometric centroid. 

Again note although it will pictorially be represented at the midpoint on the hypotenuse, 
the triangular node will actually be placed at its centroid for computations. As shown 
in Figure 19 the user responded YES, resulting in Figures 20 - 22. Note in Figure 20, 
the layer number is denoted. In this case, the user is creating a three-layer model. 
The numbering begins in the lower lefthand corner and increments to the right. The 
node centers are in line with one another and the node boundaries are identified . Note 
also the void in the model and that the software eliminated it from the sequencing. 
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FIG. 20 INTERNAL NODES - LAYER 1 
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Also note nodes 41, 42 and 48, 49. Their positions are superimposed on one 
another again only for bookkeeping. Each node will go into its respective triangle at 
the geometric centroid. Once the plot is completed, the user must hit the carriage 
return to reactivate the routine, at which time the plot will be erased. In this case, 
once the first layer is completed and erased, the second layer is plotted (Figure 21). 
Note the first node plotted is 56, one greater than the last previous node plotted in 
layer one. The nodes are sequenced exactly as previously done and the plot is labeled 
LAYER 2. Hitting the carriage return plots the third and last layer (Figure 22) , labeled 
as before, this time starting with node 111 and ending with 165. 

Refer back to the option to set the starting node number , conductance , and 
radiation links for the demonstration included in Figures 23 - 25. Here the user chose 
to start numbering the nodes at 21, the conductance block at 45, and the radiation block 
at 37. As shown in Figure 24, the first node presented in numbered 21. In the data 
set created with the actual couplings (which will be discussed fully later) , the first 
conductance number would be 45 and the first radiation number, 37. This is extremely 
useful when you are building a model in steps . 


nngen 

INPUT TOLERANCE E1E.5 

.1 

ENTER NODE N0., COND N0., RAD NO*, VOU UISH TO START UITH. IF 1 CR 
21,45,37, 

DO YOU UISH TO SEE A PLOT OF THE INTERNAL NODE CONFIGURATION*? 
yes 

2165-023B 

FIG. 23 INTERNAL NODES EXTENDED NUMBERING SCHEME 

The program then provides a rotate option for the user which allows the user to 
see graphically an orthogonal view of the network. The program queries the user as 
to how many layers he or she wishes to see and the layer numbers - (starting from 
top to bottom). As shown, the user in this case chose three layers: 1,4, and 7 
(Figure 26). The results are shown m Figures 27 - 29. Compare Figure 27 with 
Figures 20, 21, and 22. Note the rotation that took place. Similarily with Figures 
27 - 29 and the same figures. Note also that the layer numbering is from top to bottom. 
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FIG. 25 INTERNAL NODES EXTENDED NUMBERING SCHEME - LAYER 2 
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«||M 

HOU HANY LAYERS DO YOU UANT 
• 3 # 

ENTER LAYER NUNBERS 
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FIG. 26 ROTATE OPTION 



















Upon hitting a carriage return the user will erase the screen and will be asked 

whether he wants to employ the boundary option (Figure 30). hi this case, the user 

replied NO and was asked whether he wanted to create a CAVES input file. Once 

the user responds YES, the program has established the entire bookkeeping scheme 

» 

of the model in all dimensions, i.e. , it knows that node 20 (Figure 20) sees node 19 
and 21 in the x direction and node 12 and a boundary in the y direction. It also knows 
that node 20 sees node 75 (Figure 21) and a boundary in the z direction. It has com- 
puted the conductance links to all nodes (other than boundary) using the distances 
referenced from the node centers. It has computed all capacitances and established 
radiation areas (all areas to boundary). See Section 6 for a description of coupling 
computation. 

The program then asks the user for a control card and a series of title cards 
(up to four lines). If he has less than four, the user need only hit a carriage return 
for the next line to signal the program that the title entry is finished. The program 
then requests the START, DTIME, FINAL TIME, and SIGMA (Stephan Boltzman 
Constant) for execution and the initial temperature of iterated nodes. 

The program then provides the user with the ability to create the convective 

boundary links. The user will be asked successively for the convective heat transfer 

coefficient for ± x direction ± y direction and ± z direction and the heat transfer 

coefficient for any diagonal boundary surfaces (hypotenuse) . A carriage return indicates 

zero heat transfer. As demonstrated in Figure 30, the user indicated that a 100. 

2 

BTU/hr ft F coefficient was to be applied to any boundary surface facing the positive 
X direction (see Section 6). There should be a zero coupling to the negative x direction. 
For the y direction he chose a zero couplmg to positive y and 200 BTU/hr ft^ F to 
the negative y direction. The user also indicated a zero coupling in both z directions 
and any diagonal (carriage return with no response indicates zero couplings) . In 
Section 6 the user can see a description of how the couplings are created. Briefly 
each direction +x, -y, etc. , is given a specific boundary node number, i.e. , positive 
X, 993, negative x, 994; positive y, 991; etc. In checking the couplings established 
in a partial output data set (Figure 31) , note that there are no couplings to node 991 
(4-y), 994 (-X), 995 (hypotenuse), and 998 and 999 (±z). 
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DO VOU UISH TO USE THE BOUNDARV OPTION? 
no 

DO YOU UANT TO CREATE A CAUE3 INPUT FILE? 

ye5 

PKJTCD TUF fTFI 

(NUMBND,NUniTR>NUMCON,NUI1RAD^NUnNOD,NUnFLX>NSTR>NEUALU, IDIAG^ IFLUX ) 
4 , 165 , 800 ^ 243 , 999 > 0 > 15000 , 5 , 0 , 0 ^ 

ENTER TITLE CARDS 4 MAX, AFTER LAST CARD HIT CAR, RET, 
test 

of 

nngen 

ENTER START TinE,DTIME,FINAL TIME,SIGI1A 
0 ,,* 1 , 1 *, 1 , 71 40 *~ 8 , 

ENTER INITIAL TEMPERATURE OF ITERATED NODES 
50,, 

ENTER CONUECTING HEAT TRANSFER COEF FOR X DIRECTION, PLUS, MINUS 

100 , 0 , 

ENTER CONUECTING HEAT TRANSFER COEF FOR V DIRECTION,PLUS,NINUS 
0 # 200 # 

ENTER CONUECTING HEAT TRANSFER COEF FOR Z DIRECTION,PLUS,HINUS 
ENTER CONUECTING HEAT TRANSFER COEF FOR HYPOTENUSE 

2165-030B 


FIG. 30 CREATION OF INPUT FILE TO CAVE3 
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FIG. 31 SAMPLE OUTPUT DISPLAYING BOUNDARY NODES (SHEET 1 OF 2) 
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FIG. 31 SAMPLE OUTPUT DISPLAYING BOUNDARY NODES (SHEET 2 OF 2) 




Following the conduction block, the program supplies a radiation block which 
contains the link bookkeeping and the value of the node surface area exposed to the 
boundary. This block can be replaced with the actual radiation couplings when 
computed. A dummy block is supplied for Heat Sources where the user can supply 
the needed information. This block is followed by a Conductance Area Constant Block 
representing the A/ Ax term in computing the conductance link. This term is 
extremely useful when the user has a condition of changmg properties . The data 
supplied is the conductance number, the constant values from the first and second 
nodes listed for that particular conductance number, and the sector number of each 
node. An example of the output is shown in Figure 32. 

The final block of data is a dummy table to allow the user to supply one if 
needed. At this point the program will establish an output file on the system, and 
it will be in a format conductive to running CAVES. The nodal arrangement will be 
that for all mternal nodes. A segmented sample (because of its large size) is 
presented in Figure 33. 

Figure 34 represents the second option provided, i.e. , the boundary option. 

Had the user responded YES to the boundary option question, the program would have 
established all thermal parameters based on boundary node center locations. Figures 
35 through 37 represent plots of the boundary option used on the model shown in Figure 
3. A description of the boundary node conductance and orientation is described in 
Section 6. Figure 35 is identical to Figure 20 except any node exposed to a boundary 
is moved to the surface. The bookkeeping is identical but, when computing the 
conductance, a new conductive distance is used based upon its location relative to the 
node it sees (see Section 6) . The program performs identically as before and 
establishes an output data file ready for input to the CAVE3 code. 
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FIG. 32 SAMPLE OUTPUT DISPLAYING CONDUCTION CONSTANT 
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FIG. 33 SAMPLE OUTPUT FILE (SHEET 2 OF 3) 
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DO YOU WISH TO USE THE BOUNDARY OPTION'? 
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DO YOU UISH TO SEE A PLOT OF THE BOUNDARY NODE CONFIGURATION? 
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FIG. 34 BOUNDARY NODE OPTION 
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FIG. 36 BOUNDARY NODE OPTION - LAYER 2 
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Section 6 


DETERMINATION OF COUPLINGS 

The units of length and area are based on the model dimensions in the input 
file. Consistency of units is the responsibility of the user. Both conductance 
couplings and radiating areas are calculated by the program. 

6.1 CONDUCTANCE (K) 

6.1.1 Rectangular Nodes 

6. 1.1.1 Internal - The determination of conductive couplings is calculated as K = kA/L 
where A is the cross-sectional area between nodes , and L is the distance between 
centroids. 

6. 1.1.2 Boundary Option - The conductive coupling calculation is the same as that 
for internal nodes except that the additional length out to the boundary of surface nodes 
is taken into account. This approach is applied in all three dimensions. 

6.1.2 Triangular Nodes 

6. 1.2.1 Internal - The node center is assumed to be at the centroid of the triangle. 

The determination of conductance is illustrated in Figure 38. 



FIG. 38 DETERMINATION OF CONDUCTANCE LINKS - TRIANGLE TO RECTANGLE 
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For triangular nodes, the cross-sectional area is the average of the area at the 
node surface and that at the centroid parallel to the node surface . 


A 


avg 


A +A 
c s 


2 


= 5/6 Ag 


^2 


kA k, (5/6 (AY) t) ^ AYt 

= -i = 5/2 k^ 

L AXj^/3 ^ AX 

kA kg AYt 2kg AYt 

IT'AXg/s a5c 


s = surface 
c = centroid 


KiV 

%2 -K. 




Another internal triangular node configuration is illustrated in Figure 39. In 
this case the average cross-sectional area between the centroid and the sloped boundary 
is: 

A =5/6 (>/(AX)2 + (AY)2 )t 

and 

H-interface = >/(AX)^ +(AY)^ )= 1/6 >/(AX)^ +(AY)^ 

thus 

k^A k(5/6>/“ ) t 

K, 5 kt 

^ L l/6>/ 

and as before 

KiKg 

%2=iqTiq 


♦Note that the displacement between nodes in the y-direction has been neglected. 
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2165*042B 

FIG. 39 DETERMINATION OF CONDUCTANCE LINKS - TRIANGLE TO TRIANGLE 


It should be noted that this representation is accurate for isosceles triangular nodes . 
For triangular nodes with a large aspect ratio (e. g. , greater than 2:1) this rep- 
resentation may not give adequate results so that a more nearly square grid in 
desireable when using triai^ular nodes. 

6. 1.2.2 Boundary Option - the conductive coupling calculation is the same as for 
mternal nodes except that the additional length out to the surface or comer of the node 
is taken into account, hi each case the average cross-sectional area taken over the 
length L is used. Figure 40 illustrates the calculation details for the possible 
configurations . 

6.2 RADIATION 

As an aid to the user, exposed surface areas of the nodes are set up in the 
radiation couplmg block. These areas can either be modified or deleted. Seven 
fictitious boundary node designations have been used to allow rapid identification of 
node surfaces. These boundary nodes are set up as the second node number in the 
radiation coupling format. The boundary nodes are defined in the following table: 
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CONFIGURATION 


EQUATIONS 




Avg x-direction area = 

A^ = 3/4Ay *t 

t = layer thickness 

Kix = k [Aj^] = 1.5 k Ay *t 

Ax/2 Ax 

alsoK.jy = kAy = 1.5 k Ax * t 

Ay/2 Ay 


A^ = 3/4 y * t 
Ay = 3/4 X * t 

1.5 k Ay *t 
Ax 

Kiy = 1.5k Ax *t 
Ay 


A^y 2 avg= [i/2 >/(Ax)2 + (Ay)2 
L=1/2 >/(Ax)2+ (Ay)2 

*^1xy ” ^xyz avg “ 

L 
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FIG. 40 DETERMINATION OF CONDUCTANCE LINKS - BOUNDARY NODES 
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+ y side 

xz 
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- y side 

xz 
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+ X side 

yz 
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- X side 

yz 

995 

X - y (hypotenuse) 

xyz 

998 

- z side 

xy 
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Section 7 


BOUNDARY NODE ORIENTATION 


When the boundary option is employed , the internal node configuration is 
altered to a prescribed surface configuration. The actual location of the surface 
node depends on the boundary orientation with respect to the node in question. Figure 
41 displays the movement of the internal node to the surface. The figure is broken 
down into 6 (A - F) representative surface nodes. "A" represents an internal node 
surrounded on all sides but one. The node is then placed at the midpomt of the surface 
exposed to boundary when the boundary option is activated. Although only one orien- 
tation of "A" is shown, the movement could have been in any direction, i.e. , the 
boundary could have been on the bottom (in which case the node would have been placed 
there), etc. 

"B" represents a node surrounded on two adjacent sides and the node is placed 
at the comer. ”C" represents a triangular node the hypotenuse of which sees boundary. 
The node is placed at the midpoint of the hypotenuse. "D" displays a triangular node 
bases of which are exposed to boundary , in which case the node is placed at the 
comer. "E" and "F" represent triangular nodes where only one base is exposed to 
boundary, in which the node is placed at the midpoint of that base. 

Figure 42 displays improperly designed nodes that cannot employ the boundary 
option. When any of these arrangements are encountered, the program will leave 
the mternal node unaltered and write a message m the output file warning the user 
that the node was not altered. 










This Page intentionally Left Blank 



Section 8 


REFERENCES 


1. J. Palmieri and K. Rathjen, "CAVES; A General Transient Heat Transfer Computer 
Code Utilizing Eigenvectors and Eigenvalues," NASA CR145290, Feburary 1978. 


63 



2. Government Accession No 


3 Recipient's Catalog No 


1, Report No 

NASA CR-l 65860 

2. Government Accession No 

3 Recipient's Catalog No 

4 Title and Subtitle 

5 Reoort Date 

December 1 982 



6 Performing Organization Code 

7. Author(s) 

Joseph V. Palmier! 
Kenneth A. Rathjen 

8 Performing Organization Report No 

10 Work Unit No 

9 Performing Organization Name and Address 



Grumman Aerospace Corporation 
Bethpage, N.Y. 11714 

11 Contract or Grant No 

NAS 1-15367 



13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 


Contractor Report 

National Aeronautics Space Administration 
Washington, D.C. 20546 

14 Army Project No 


IS. Supplementary Notes 

Langley Technical Monitor: James L. Hunt 


16 Abstract 


Under contract NASA-14643, Grumman developed the computer code CAVES (”A 
General Transient Heat Transfer Computer Code Utilizing Eigenvectors and Eigenvalues”) 
which performed a specialized transient thermal analysis. 

Under the present funding, a new extension of CAVE3 code was developed that 
automates the creation of a finite difference math model in digital form ready for input 
to the CAVES code. The new software. Nodal Network Generator, is broken into two 
segments. One segment generates the model geometry using a Tektronix Tablet Digitizer and 
the other generates the actual finite difference model and allows for graphic verification 
using Tektronix 4014 Graphic Scope. 

This report describes how to use the Nodal Network Generator. 


17 Key Words (Suggested by Author(sM 

Aerodynamic Heating 
Transient Heat Transfer 
Eigenvalue-Eigenvector Solution 
Structural Temperature 
Hybnd Analytical Numerical Method 

18 Distribution Statement 

Unclassified — Unlimited 

19 Security Oassif (of this reporti 

UNCLASSIFIED 

20 Security Classif (of this page! 

UNCLASSIFIED 

21 No of Pages 

73 

22 Price* 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley • 1979 
























End of Document 



